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Abstract 


We study reheating in theories where inflation is completed by a first-order phase 
transition. In these scenarios, the Universe decays from its false vacuum state by bubble 
nucleation. In the first stage of reheating, vacuum energy is converted into kinetic energy 
for the bubble walls. To help understand this phase we derive a simple expression for the 
equation of state of a universe filled with expanding bubbles. Eventually, the bubble walls 
collide. We present numerical simulations of two-bubble collisions clarifying and extending 
previous work by Hawking, Moss, and Stewart. Our results indicate that wall energy is 
efficiently converted into coherent scalar waves. We go on to discuss particle production 
due to quantum effects. These effects lead to the decay of the coherent scalar waves. In 
addition, they lead to direct particle production during bubble- wall collisions. We calculate 
particle production for colliding walls in both sine-Gordon and <p 4 theories and show that 
it is far more efficient in the <f > 4 case. The relevance of our work for recently proposed 
models of first-order inflation is discussed. 
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1. Introduction 


In Guth’s original version of inflation [1], the Universe decays from its initial false- 
vacuum state by the nucleation of true-vacuum bubbles. As the bubbles expand, the 
energy in the fake vacuum is converted into kinetic energy for the bubble walls. In theory, 
the walls eventually collide, the true vacuum percolates, and wall energy is converted into 
radiation, reheating the Universe. Unfortunately, in Guth’s model the nucleation and 
expansion of bubbles cannot keep up with the exponential expansion of the regions still 
trapped in the false vacuum, and the true vacuum never percolates [2], This is known in 
the literature as the ‘graceful exit’ problem. 

Clearly, first-order inflation (generically, any model in which inflation is completed 
by a strongly first-order phase transition) can work if either the nucleation or expansion 
rate changes during the inflationary phase. For example, in the extended inflation scenario 
proposed by La and Steinhardt [3] the Hubble parameter in the vacuum-energy dominated 
( p = — p) Universe decreases with time, thus enabling the true vacuum to percolate. While 
La and Steinhardt’s proposal suffered a graceful exit problem of its own [4], there has been 
no shortage of theories which purport to save the general idea of first-order inflation [5]. 

Given the renewed interest in first-order inflation, we felt it an opportune time to 
examine some of the nuts and bolts of reheating in these models. Reheating in new [6] and 
chaotic [7] inflation has been examined in some detail [8]. The analysis in these scenarios 
is simplified by the fact that the Universe can be treated as homogeneous (Friedmann- 
Robertson- Walker) throughout. Inflation occurs when some order parameter or scalar 
field finds itself displaced from the minimum of its potential. At first, the field slowly 
rolls towards this minimum. During this phase the energy density p of the Universe is 
dominated by vacuum energy with p = const' = p„ ac ; the pressure p = -p; and the scale 
factor a oc e H,t , where Hj = 87rp vac /3m 2 l . (We use units where h = c = kg = l and 

= G N ^t on ~ 10 19 GeV is the Planck mass). As the inflaton reaches the minimum, 
the slope of the potential increases and the field begins to oscillate about the minimum 
on a time scale short compared to the Hubble time. These oscillations behave at first like 
pressureless nonrelativistic matter so that a oc t 2 ^ 3 . Eventually, the oscillations decay into 
relativistic particles and the Universe becomes radiation-dominated. We note that if the 
lifetime of the inflaton is long compared to HJ 1 , most of the energy in the <f> field will be 
diluted by cosmological expansion and the reheat temperature Trh ( ~ (mp/T^) 1 / 2 ) will 
be well below piac, the energy scale associated with the false vacuum. 

The physics of reheating in first-order inflation is essentially the same as in new and 
chaotic inflation: Energy initially stored in a coherent scalar field must be converted into 
radiation. However, the situation is considerably more difficult to analyze. Here, both 
the field <f> and the metric are inhomogeneous. Furthermore, reheating involves a mix of 
quantum physics (bubble nucleation and particle creation) and classical physics (expansion 
and collisions of bubbles). Difficulties aside, a number of authors have focused on the 
potentially rich phenomenology associated with these scenarios. For example, it has been 
suggested that gravitational waves [9], black holes [10] [11], topological defects [12], and the 
baryon asymmetry [13] [11] may have been produced during the phase transition. Whether 
or not such phenomena actually occur depends in part on the details of reheating. For 
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though any asymmetric double-well potential will do. We shall generically refer to models 
of this type as <£ 4 models. Focus for the moment on the case where e< Aso that the energy 
difference between true and false minima is small compared to the height of the barrier. To 
lowest order in e, the true and false minima are at — <p 0 and <f > 0 respectively. The mass of 
^-excitations in either of these minima is given by m $ = A 1 / 2 <t> 0 , and p vac = 2 e<p*. e < A 
corresponds to the thin-wall limit in which the radius of the bubbles is much greater than 
the thickness of the wall separating the interior true vacuum region from the exterior false 
vacuum region. One nice feature of the thin-wall limit is that the properties of individual 
bubbles can be determined analytically. In particular, the thickness of the bubble wall is 
A = 2/ (A^Vo) = 2 m^ 1 , the surface energy density is o — 2A 1 / 2 </> 3 /3, and the radius of 

the bubble at the time of nucleation is R 0 = A 1 / 2 / (e(f> 0 ) = 3 cr/p vac - Ro/& — A/e 1 so 
the bubble walls are indeed thin. 

The results quoted above are valid so long as gravity is not important. Gravity be- 
comes important when the radius of the bubbles becomes comparable to the Hubble radius 

[16]. For the case at hand, the Hubble radius is given by H~ l = \S>T:p V ac/^ m pi I — 

m p i / (e 1 / 2 ^ 2 ). RH ~ (A/e) 1 / 2 (t> 0 /m p i is therefore a measure of how important gravity is. 

The nucleation rate per unit volume (number of bubbles nucleated per unit four- 
volume) is [17] 

r ~ A (m^) 4 e _5e (2.3) 

where A is a constant of order unity and Se is the Euclidean action for the bounce solution 
corresponding to a critical bubble. For the thin wall case, Se = 7r 2 A 2 /3e 3 . Perhaps the 
most important parameter describing a first-order inflation theory is 77 = TJH 4 , the so- 
called percolation parameter. 77 gives roughly the number of bubbles nucleated per horizon 
volume per Hubble time. For the case at hand 



The requirement that inflation last long enough to solve the horizon and flatness problems 
constrains 77 to be less than about 10 -3 , whereas percolation of the true vacuum bubbles 
can only occur if 77 > 0(0.1). In Guth’s model, 77 is constant and the model is untenable. 
In the models of Ref. [5] , 77 varies in time; a period of inflation with 77 < 10 -3 precedes the 
percolation phase in which 77 is large. 

The mean separation between bubbles is roughly D — T -1 / 4 = 77 -1 / 4 # -1 ~ 
(A 1 / 2 ^) 1 e^ E / 4 . Certainly, a key unknown is the value (or range of values) 77 has when 
percolation occurs. (Remember that 77 is increasing with time.) For 77 ^ 1, D will be less 
than the Hubble radius and there will be many bubbles within a given Hubble volume. In 
this case, cosmological expansion can be neglected in treating the dynamics of the bubbles. 
On the other hand, if 77 ~ 1 when percolation occurs, D will be of order the horizon and 
cosmological expansion will play an essential role in the bubble dynamics. In what follows 
we will neglect, for the most part, cosmological expansion in treating the bubbles, but only 
because it simplifies the calculations. We leave the full problem of bubble dynamics in an 
expanding spacetime for future investigations. 
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Let TZ be the stress energy for the wall. In the thin wall approximation TZ has a ,5- 
function singularity across the wall. It is convenient to define the quantity 

S o6 = Jdl TZ (3.3) 

where l is the proper distance through the surface in the direction of £ a . We note that 
dl = 7 dr. 

The three metric intrinsic to the wall is 

h-ab — 9ab £a£b • (3-4) 

For domain walls, S„6 = —ah a i,. In addition bubbles nucleated by quantum tunneling have 
3 a/R 0 = p vac - It is straightfoward to show that 

Q -1 J d z xTZ 


ft -1 J d z xTZ 

Combining Eq.(3.2) and Eq.(3.5) we find 

P — Pvac 
P ~ Pvac fl 

For R < D, p ~ -p vac as expected; the Universe in this case is essentially dominated by 
vacuum energy. If, on the other hand, bubbles are nucleated with R 0 ~ D, then p = -2p/3 
right after nucleation. This is just the equation of state for a wall-dominated Universe. 
In fact, the Universe is still undergoing power-law inflation when the bubbles first collide 
[20]. Finally, if R 0 <C D, the walls will be relativistic by the time they collide, and the 
equation of state just before bubble collisions will be p = p/3. 


ZjaR 2 

D 2 

R z 

~ Pvac~£$ ( 3 * 5 ) 

“ ^ ac £)3 \^ 7 2 3 J ■ 
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4. Bubble Collisions: Classical Treatment 


Bubble collisions provide the next step in reheating. Collisions release energy bound in 
the walls through both quantum and classical processes. Quantum effects will be discussed 
in Section 5. Here, we treat the bubble walls as classical field configurations and show that 
classical scalar waves are emitted during a collision. We explore this process by studying 
the collision of two expanding bubbles. 
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For |r| < |f| (Region II), we require a different coordinate system {p\ ip' ,9,p>) with 

x — p 1 sinh ip' sin 9 cos ip z = p' sinh ip' cos 9 

y = p' sinh ip' sin 9 sin ip t = p' cosh ip' . 

In order to find the field configuration in this region one needs to take the field configuration 
on the p = 0 hypersurface (from the Region I solution) and evolve into Region II. Since 
the field configuration only depends on p this amounts to solving an ordinary differential 
equation. For a single thin-wall bubble, the bubble wall is entirely in Region I. The field 
is essentially constant in Region II, so that evolving the field is usually unnecessary. 

The situation changes with more than one bubble. Consider the nucleation of two 
bubbles. In general this is a complicated process and one that has yet to be treated in 
the literature. If, however, the bubbles are widely separated at the time of nucleation, 
then they can be treated as noninteracting (the dilute instanton approximation), and the 
generalization from the single bounce solution is straightfoward. For two bubbles, the axis 
joining their centers will be a preferred direction, so that the solution to the Euclidean 
equations of motion for noninteracting bubbles will have 0(3) symmetry. We consider 
first a coordinate system in which the r = 0 hypersurface intersects the centers of the 
two bubbles. The field configuration on this hypersurface becomes initial data for two 
simultaneously nucleated expanding bubbles in Minkowski space, and one could solve for 
the field configurations numerically by solving the classical equations of motion. As in the 
single bubble case, a more efficient technique for finding the evolution exploits the extra 
boost symmetry of the problem. Here the 0(3) symmetry of the two-bubble Euclidean 
bounce translates to an 0(2, 1) symmetry in Minkowski space. Let the 2 -axis correspond 
to the the line connecti ng the ce nters of the bubbles. As before, spacetime is divided into 
two regions. For |t| < yfx* + y 2 we choose the coordinates ( s,ip,9,z ) with 

x = s cosh ip sin 9 z = z 

(4 3) 

y = s cosh ip cos 9 t = s sinh ip . v ' 

The solution in this region of space time is g iven by the analytic continuation of the two 
bubble bounce solution. For |t| > \Jx 2 + y 2 we take 

x = s' sinh ip' sin 9 z = z 

y = s' sinh ip' cos 9 t = s' cosh ip' . ^ ^ 

To find the field configuration in this region, we take the solution on the t = r = v /^ 2 TjT ^2 
hypersurface and solve the equations of motion. Here the equations we are required to 
solve are 1 + 1 (s and z) partial differential equations. For the two bubble case, all of the 
interaction between the bubbles takes place in this second region. 

The procedure outlined above was developed and used by HMS. Here we consider 
bubbles in a theory described by Eqs.(2.1,2.2) with e = 0.1. In Fig. la, we show the field 
<t>(s,z) in the collision region for two simultaneously nucleated bubbles. In Fig. lb, the 
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an observer at z = 0 but arbitrary r, r(r) = yjt 2 f + r 2 - \/tj + r 2 . For r < {t x ,t f ), 

T ( r ) = To (l - r 2 /2tit f ) whereas for r » (<*,</), r(r) = r 0 (f* + t f ) /2 r. In either case, the 
time interval decreases with increasing r, a result which can be read off of Fig. 2. 

We have carried out a variety of numerical studies and found that generally the radia- 
tion of scalar waves is reasonably efficient. For example, in the simulations described above 
we find that most of the energy in the center of the two-bubble system is radiated away 
after only a few collisions. This result holds for a wide range of initial separations. We 
have also simulated the collision of two infinite plane walls in a theory with a symmetric 
double well potential (Eq.(2.2) with e = 0) and find that the percentage of energy radiated 
during a single collision is roughly constant and ~ 30-40% for a wide range (0.3 < 7 < 10) 
of initial velocities (see, for example, Ref.[21]). 

One can also consider first-order inflation in sine-Gordon (SG)-type theories. Here 
the potential is of the form 

V = (cos + VM)) (4.5) 

where there are N local minima at <j> = 2 n<p 0 and V A creates a small asymmetry among 
these minima. For the most part, the dynamics of wall- wall interactions is determined 
by the first term in V. As is well known, SG kinks in a theory with V A = 0 are true 
solitons, and infinite plane-symmetric walls pass through one another without dissipating 
any energy. However, if the walls are curved (or if V A ± 0) they will produce scalar 
radiation [22], though not as efficiently as <p 4 walls. Reheating, in fact may be very different 
in these models (see also Section 5), though no detailed work has yet been done. 


5. Particle Production 


In order to reheat the Universe, the inflaton must couple to ordinary particles. This 
will allow the classical scalar waves described above to decay, eventually filling the Universe 
with a thermal bath at a temperature Trh. These couplings also lead to direct production 
of particles during collisions between walls. 

In this section we discuss this direct particle production. In particular, we con- 
sider production of fermions arising from interactions of the form C\ = g\<pipip and 
£2 = 92f- 1 i>'Yvi' t Pd* t 4>, where / has units of mass. C\ is the typical Yukawa coupling of a 
scalar field to fermions. £ 2 -type couplings arise if <j> is a Goldstone (or pseudo-Goldstone) 
boson. In this case, the potential for <p is of the sine-Gordon (SG) type. As we will see, 
particle production in these two theories is dramatically different, though this is due as 
much to the peculiar properties of the SG walls as to the different form of the coupling 
between <j> and ip. 

Our philosophy is to treat </> (the bubbles or walls) as a classical, external field and the 
fermions as quantum fields in the presence of this source (see, for example, Ref. [23]). In 
so doing we make no attempt to treat backreaction of particle production on the evolution 
of the walls. 
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As a first example, consider an infinite plane-symmetric domain wall. To keep things 
simple we assume a model with degenerate vacua so that a noninteracting wall will move 
with constant velocity. Since the result must be Lorentz-invariant, we can work in the rest 
frame of the wall. Clearly, in this frame po = 0 for all modes, so that p 2 < 0 and there is 
no particle production. 

Next, consider the head-on collision of two plane-symmetric walls. The walls can come 
from either <f > 4 or SG theories. As we will see, the crucial difference between the two models 
is that <j > 4 walls scatter off one another while SG walls pass through one another. 

In principle, the calculation should be straightforward. One starts with the field 
configuration describing a collision. In the case of plane-symmetric walls moving in the 
z-direction, <f>(x,t) = <fr(z,t). The next step is to calculate the Fourier transform. Here 
= ( 2Tr) 2 6(k x )6(k y )<j>(k z ,u) ). The final step is to substitute into Eq.(5.7) and inte- 
grate over p. This will give the number of particles produced per unit area: 

J = 2 J Im (r (2) (^ 2 - * 2 )) (5.8) 


where k = k z . 

The field configuration for colliding </> 4 -walls can be found numerically. However, the 
result will be quite complicated, as the scattering of two 4 > 4 walls is inelastic. In general, 
scalar waves are produced in the collision whose decay will contribute to the particle 
production. It is difficult to untangle this contribution from that of the actual collision. 
We therefore choose to model the collision by treating the walls as infinitesimally thin and 
assuming that the scattering is perfectly elastic. This model should give accurate results 
for k < 7771^ and u 'ym^/v in the case of relativistic walls. 

Our ansatz for the field configuration is 


( 00, 

for vt < z < —vt and t < 0; 

0 

II 

for — vt < z < vt and t > 0; 

v -<t> 0, 

otherwise 


and the Fourier transform of the field is 


4>(k,u) 


8u0 o 

OJ 2 — k 2 v 2 


(5.10) 


The slow power law-fall off at large u> 2 - k 2 v 2 is due to the fact that in our ansatz, $ 
is discontinuous. For realistic “thick” domain walls, 4> would cut off exponentially for 
k > 77M0 and u > Substituting Eqs.(5.5,5.10) into Eq.(5.7) we estimate 


N 

A 



(5.11) 


where we have assumed that 7777^ > 2p and 7 > 1. The energy per unit area radiated by 
the walls is 


E 

A 




(5-12) 
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The final result can be written in the form 


4>(u,k) 




Wo Sinh (2^) 

cosh (r ~g ) - cosh 

V ) \l rn * ) 


(5.20) 


This result is similar in many respects to Eq.(5.10). For u < 7 vm# and k < 7 m*, 

<f> oc (a; 2 — v 2 k 2 ) .In addition 4> is exponentially damped for larger values of k and u>, as 
anticipated above. We estimate that the number of particles produced per unit area is 



and the energy lost per unit area is 



(5.21) 


(5.22) 


where again we have assumed that ym^ 2fi. This represents a fraction (<?2 m 4>/7/) 2 
of the total energy in the walls and is less than that found in the <f > 4 case by a factor 

i m 4>/ f) 7 • The (m^//) 2 term comes from the nature of the coupling C 2 and is usuallv 

quite small. 

Particle production in this case is due to a slight slowing down of the walls during 
their interaction. As 7 becomes large, the forces between the walls are able to effect a 
smaller change in velocity, leading to a smaller result for N/A. In principle, there will be 
similar effects in the 4 > 4 case leading to corrections of order l/ 7 n to Eqs. (5.11) and (5.12). 

Evidently, we can produce particles up to energy 7 m. If the bubble walls are highly 
relativistic when they collide, there will be the possibility of producing particles well above 
the mass of the infiaton. 


6. Conclusions 


It is relatively easy to find a cosmological scenario which inflates. The challenge is 
in bringing an end to inflation so that the resultant universe resembles the one we live 
in. The vacuum energy which drives inflation must be converted into relativistic particles. 
Furthermore, fluctuations which arise during inflation should not lead to unacceptably 
large distortions in the microwave background. It is a bonus if they provide the seed 
perturbations necessary to drive the formation of large-scale structure. In addition, relics 
of the inflationary epoch such as gravity waves and topological defects may have survived 
until the present epoch. If detected these phenomena would provide a window to the very 
early Unverse. 

The search for a truly graceful exit from inflation has had limited success. As is 
well known, most new and chaotic inflation scenarios require ‘fine tuning’ to satisfy the 
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black, the true vacuum in light gray. Regions where <f> < -<f> 0 appear white. In Fig. 
lb the initial separation beteen the bubbles has been doubled and the coordinates 
have been scaled by a factor 1 /2. 

2. Spacetime diagram in the t - r plane for the two-bubble collision shown in Fig. la. 
For this diagram, 2 = 0. 

3. Field configuration shown at different times during the collision. Note the scalar 
radiation emanating from the collision region. 

4. Spacetime diagrams illustrating a big bubble hitting a small bubble. Fig. 4a shows 
the two-bubble collision in the equal bubble (simultaneous-nucleation) frame Here 
the radiation is symmetric about 2 = 0. Fig. 4b shows the same configuration in a 
Lorentz-boosted frame. The bubble on the left is nucleated first and is larger at the 
time of the collision. In this frame, all of the radiation is moving to the right, away 
from the nucleation site of the big bubble. 
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